Measurements of Quasiparticle Tunneling Dynamics in a Bandgap-Engineered Transmon Qubit 
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We have engineered the bandgap profile of transmon qubits by combining oxygen-doped Al for tunnel junc- 
tion electrodes and clean Al as quasiparticle traps to investigate energy relaxation due to quasiparticle tunneling. 
The relaxation time T\ of the qubits is shown to be insensitive to this bandgap engineering. Operating at rela- 
tively low EjlE(j makes the transmon transition frequency distinctly dependent on the charge parity, allowing 
us to detect the quasiparticles tunneling across the qubit junction. Quasiparticle kinetics have been studied by 
monitoring the frequency switching due to even/odd parity change in real time. It shows the switching time is 
faster than 10 lis, indicating quasiparticle-induced relaxation has to be reduced to achieve T\ much longer than 

100 ;US. 



Quantum information processing based on superconduct- 
ing qubits has made tremendous progress towards realizing a 
practical quantum computer in the last few years . How- 
ever, the coherence times of superconducting qubits still need 
to improve to reach the fault-tolerant threshold. For exam- 
ple, the single-qubit gate error rate is limited by qubit deco- 
herence Jj-kjfl. Understanding decoherence mechanisms, in 
particular those responsible for qubit relaxation, is therefore 
crucial. Quasiparticles have received significant attention re- 
cently as one such possible limiting factor f7H9[]. At low tem- 
peratures, thermal-equilibrium quasiparticles should be irrel- 
evant because their density is exponentially suppressed. In 
practice, however, non-equilibrium quasiparticles are present 
from unknown sources 111 (A ll 111 . The key question is then 
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whether these non-equilibrium quasiparticles are currently 
limiting qubit relaxation, and if not, what limit they will ul- 
timately impose. The answer will be relevant to all supercon- 
ducting qubits. 

Quasiparticles in qubit electrodes do not themselves cause 
significant qubit decoherence. Instead, it is the dissipative 
and incoherent tunneling of quasiparticles across a Joseph- 
son junction that leads to decoherence. The tunneling pro- 
cess is characterized by the quasiparticle current spectral den- 
sity S qp (co). The low frequency component in S qp (co) can 
cause dephasing, but this dephasing channel can be elimi- 
nated b y re ducing the qubit's sensitivity to charge noise, as 
in flux 113, [l3ll . phase 113], and transmon lfl5ll qubits. The 
high frequency component will cause energy relaxation in all 
types of superconducting qubits. In a transmon qubit, which 
has the energy level structure shown in Fig. [Ha), quasiparticle 
tunneling leads to two observable rates: the qubit relaxation 
2 5„„(cooi ) and the odd/even parity switching 
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rate = \M oe \ 2 S qp (co oe ), where corji /27T and co oe /2n are the 
qubit ground-to-excited and odd-to-even transition frequen- 
cies respectively, and Mq\ and M oe are the matrix elements 
describing the interaction between the qubit and quasiparti- 
cles Jlrjll . 
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FIG. 1: (a) Low-energy level structure of a transmon qubit. Sq and 
8\ : energy difference between odd/even parities for the ground and 
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is the charge dispersion. 1%% is the odd/even parity switching rate. 
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coqi S> 5i. (b) SEM image of a typical device fabricated by 
three-angle evaporation to engineer the bandgap profile, (c) Expected 
bandgap profile near the tunnel junction. 



Several approaches may be taken to reduce the detrimen- 
tal effects of quasiparticles on qubits. One obvious way is to 
reduce the density of quasiparticles by lowering the temper- 
ature, increasing the superconducting bandgap, or suppress- 
ing the generation of quasiparticles. Even if the quasiparticles 
cannot be eliminated, S qp could still be reduced by prevent- 
ing tunneling of quasiparticles. For example, the bandgap 
profile could be engineered so that quasiparticle tunneling 
becomes energetically suppressed. This bandgap engineer- 
ing approach has been successfully applied to reduce r^ 5 in 



single-Cooper-pair transistors (SCPTs) ll 1 OL 17 - 1^1 and super- 
conducting charge qubits [ 20] . 

In this paper, we adopt bandgap engineering to study the 
relaxation of transmon qubits due to quasiparticles. We mea- 
sure the qubit relaxation time T\ as a function of temper- 
ature of transmon qubits which have significantly different 
bandgaps. At low temperatures, the saturation of T\ of both 
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transmons at approximately the same level suggests a mech- 
anism other than thermal quasiparticles for qubit relaxation. 
To investigate if non-equilibrium quasiparticles are limiting 
T\, we then study the quasiparticle kinetics in a bandgap- 
engineered transmon operated in the low EjIEq regime, where 
Ej and Eq are the Josephson and the charging energy, respec- 
tively. Qubit spectroscopy shows two qubit transition frequen- 
cies associated with the even- and odd-charge states, demon- 
strating the presence and tunneling of quasiparticles across 
the qubit junction. This indicates that with our design, non- 
equilibrium quasiparticles have not been removed by bandgap 
en g inee ring , contrary to the results of experiments in SCPTs 
lild.ll8lll9ll . We study the parity switching rate T^f in the time 
domain finding 1 /T$e < 10 jus. For typical device parameters, 
the expected ratio of parity switching rate to quasiparticle- 
induced qubit relaxation rate, TfelYf^ - 10 - 100 {u&. 
Therefore, while we cannot establish quasiparticle tunneling 
as the dominant source of energy relaxation in our devices, our 
bound indicates that reducing the quasiparticle-induced decay 
rate will be necessary to achieve T\ much longer than 100 jus. 

Our transmon qubits are measured using a coplanar waveg- 
uide cavity in a conventional circuit quantum electrodynamics 
architecture [21]. All devices are measured in a cryogen-free 
dilution refrigerator with a 20 mK base temperature. Oxygen- 
doped Al, which has an energy gap A « 280 jueV (T c wl.9K) 
about 60% higher than clean Al (A w 180 jUeV, T c » 1.2 K), 
is used as the electrodes of Josephson tunnel junctions de- 
posited by standard double-angle evaporation 02211 . The oxy- 
gen dopants are introduced with a continuous O2 flow during 
the Al deposition. The same technique has also been used in 
SCPTs to realize a large bandgap [10, 18]. To create quasipar- 
ticle traps, a third layer of clean Al is deposited to cover the 
whole oxygen-doped Al layers to within ^100 nm from the 
junctions [Fig. [Tib)]. Figure[TJc) shows the expected bandgap 
profile near the junction. This profile is expected to trap quasi- 
particles away from the tunnel junction. The energy gaps of 
oxygen-doped layers are determined by independent T c mea- 
surements (not shown) of thin films evaporated under the same 
nominal conditions as for the tunnel junction. 

However, our bandgap engineering does not appear to affect 
T\ (see Table I in Supplemental Material for list of devices). A 
comparison of T\ as a function of temperature for two repre- 
sentative devices is shown in Fig.|2a): one qubit is fabricated 
with clean Al only and the other with oxygen-doped Al with- 
out the third layer quasiparticle trap. The decrease of 7\ with 
temperatures above 150 mK for the red and 250 mK for the 
blue curve, respectively, indicates the effect of thermally gen- 
erated quasiparticles I3II6J]. The higher corner temperature 
confirms that oxygen-doped Al indeed has a larger energy gap. 
The saturation of T\ at low temperatures for both clean and 
oxygen-doped Al devices at the same value indicates that ther- 
mal quasiparticles are not limiting T\. Figure [2b) shows the 
saturated T\ as a function of frequency for four qubits, each 
of which has its own flux bias line allowing individual tuning 
of the qubit frequency [2], fabricated with both oxygen-doped 
Al and quasiparticle traps described earlier. The qubit life- 
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FIG. 2: (a) Measured qubit relaxation time T\ as a function of tem- 
perature for two transmon qubits: one qubit is fabricated with clean 
Al only with a transition frequency /01 = 4.25 GHz and the other 
with oxygen-doped Al but without the third layer of quasiparticle 
traps with /01 = 5.16 GHz. Dashed black and magenta lines are 
theory for relaxation due to thermal quasiparticles, assuming A ~ 
180 fleV and 280 /ieV, respectively |8]. The solid green and cyan 
lines represent the sum of the theoretical expectations for thermal- 
equilibrium quasiparticles and the best-fit saturated T\, (b) 7) vs 
frequency on qubits fabricated with oxygen-doped Al and quasipar- 
ticle traps. Dashed red line: Purcell-induced relaxation time; solid 
blue line: Q = 65, 000; solid black line: a combination of the Purcell 
effect and a constant Q. Neither oxygen-doped Al nor quasiparticle 
traps improve T\ qualitatively. The black arrow indicates the cavity 
frequency. 
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FIG. 3: (a) Spectroscopy of a qubit with engineered bandgap as a 
function of gate-induced charge n g . The gate voltage is applied to 
the qubit through a bias tee at the cavity input port. For each pixel, a 
gaussian pulse (a = 20 ns, corresponding to a % pulse on resonance) 
is applied at the indicated frequency and the qubit is immediately 
measured. Each pixel is average of 5000 repetitions (50 ms). Darker 
pixels correspond to higher homodyne readout voltages which are 
proportional to the probability of the qubit in the excited state. An 
"eye"-shaped pattern indicates charge-e jumps associated with the 
tunneling of non-equilibrium quasiparticles. (b) Cross-section of (a) 
at rig — 0.5 [black arrow in (a)]. The charge dispersion is 60 MHz. 
We refer to the lower (upper) frequency branch as the even (odd) 
parity peak. 



times are limited to a quality factor Q ~ 65,000, similar to 
previously reported T\ on qubits fabricated with clean Al only 

Non-equilibrium quasiparticles have been observed in the 
bandgap-engineered devices as will be shown later, so here we 
will only focus on the effect on T\ from those non-equilibrium 
quasiparticles. The lack of qualitative improvements of T\ 
suggests two possible scenarios: either T\ is not limited by 
non-equilibrium quasiparticles, or the quasiparticle contribu- 
tion to T\ has not been affected by our bandgap engineering. 
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In either case, we wish to find the quasiparticle-induced qubit 
relaxation rate rf^, . However, this rate cannot be measured 
directly in the presence of other sources of relaxation, so we 
instead measure the dynamics of quasiparticle tunneling 
in the time domain, from which we can estimate Vf^. 11611 . 
In order to measure X% , we operate the transmon qubit in the 
low Ej/Eq regime where the qubit spectrum will have two 
distinct parities caused by quasiparticle tunneling across the 
junction. Thus, Y q ( ,e can be studied by monitoring the dynam- 
ics of one particular parity in real time. 

We have fabricated a single transmon qubit with both 
oxygen-doped Al and quasiparticle traps, and operate it at 
Ej/Ec ~ 12.5 by tuning the qubit frequency to /oi « 4.5 GHz, 
where Ti = 2 /is. To protect the qubit from spontaneous emis- 
sion J2J through the low- Q cavity (Q = 500), a Purcell filter 
is also integrated on chip [24]. At this large detuning from the 
cavity frequency f c = 8.072 GHz, a direct readout of the qubit 
state is difficult because of the weak dispersive interaction be- 
tween the qubit and the cavity. Instead, to enhance readout 
fidelity, the qubit is pulsed to /bi ~ 7 GHz prior to measure- 
ment. By making use of the high-power Jaynes-Cummings 
readout 02511 . we achieve a single-shot fidelity F — 42%. 

The presence and tunneling of non-equilibrium quasiparti- 
cles is demonstrated by the clearly observed "eye" pattern in 
the qubit spectroscopy as a function of gate-induced charge n g 
[Fig- 13a)] 111]]. The two peaks in Fig. 0b) have almost equal 
height, implying no preferred parity, and both parity switch- 
ing times x (odd-to-even) and % e (even-to-odd) shorter than 
the 50 ms averaging time at each pixel. 

To study the dynamics of quasiparticle tunneling, a selec- 
tive 7t pulse is applied repeatedly every 10 /is at one of two 
parity frequencies (for simplicity, we choose the odd parity 
frequency for all data presented here), immediately followed 
by a measurement of the qubit state [see Fig.|4ja)]. This selec- 
tive n pulse will excite the qubit only when the parity is odd. 
We note that due to the smallness of 8q [Fig.QJa)], the selec- 
tive n pulse will also cover the transition |0,eve«) to \\,odd), 
but for non-degenerate quasiparticles the probability of this 
unwanted transition is small compared to that of the direct 
photon transition \0,odd) to \l,odd). Thus, by interrogating 
the qubit state after the n pulse, the qubit parity can be in- 
ferred. To ensure that the qubit begins in the ground state and 
that the qubit frequency has stabilized after the rapid tuning 
from 7 to 4.5 GHz at the beginning of each experimental cycle, 
the repetition rate is set to t„ = 10 /is = 5T\. Because quasi- 
particles tunnel across the qubit junction randomly, a random 
telegraph signal (RTS) is expected. To overcome low single- 
shot readout fidelity, we perform a Fourier transform of the 
time domain data to study the power spectral density (PSD) to 
better extract the RTS information. Background charge mo- 
tion limits the experiment duration since the qubit frequency 
shifts noticeably every few minutes. In total, twelve million 
measurements are recorded and Fourier transformed. 

The stability of the readout is particularly important for this 
measurement, and so to minimize the drift during the experi- 
ment, the readout result is digitized by thresholding the mea- 
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FIG. 4: (a) Schematic of the measurement. A selective 71 pulse is 
applied to one of the two parity peaks in Fig.[3]followed by an imme- 
diate readout of the qubit state at about 7 GHz where readout fidelity 
is improved. Lines: the flux pulse sequence. The process is repeated 
every 10 /is. (b) Histogram of the readout at the odd parity peak in 
Fig. [3jb). A threshold V,/, = 19 mV is chosen to digitize the read- 
out. Inset: Schematic of an imperfect readout with false positives 
(negatives) £o (£i). 

surement results [Fig.HJb)]. From each measurement, one bit 
of information is extracted. False positives £o and negatives 
£i of the readout [Fig. |4|b) inset] will reduce the RTS am- 
plitude. Note that due to the switching of the qubit between 
the two parities, £o and £\ cannot be extracted directly from 
Fig-ISb), but can easily be obtained by combining histograms 
of the readout at the even parity peak in Fig. [2b) (see Supple- 
mental Material for details). We find £o ~ £\ = 0.29, as well 
as the probability of the qubit at the odd parity = 56%. 
Thus the readout fidelity is F = 1 - £ - E\ « 42%. 

We first test the sensitivity of the experiment to fluctuations 
of parity by applying a "% mask" to the measurement system 
to imitate the success and failure of a % pulse on the qubit. 
The n mask, generated by a field programmable gate array 
(FPGA), is a pseudo-random control sequence which enables 
or disables a % pulse applied to the qubit immediately prior 
to readout. It generates an RTS with a specified time constant 
X K = t k q + %i and a 50% duty cycle {t k q = T^i). If the parity 
switches fast compared to x K (as will be confirmed later), the 
RTS amplitude in the readout is reduced to A c0 „ — P oc uiF due 
to the parity duty cycle and the finite readout fidelity [Fig. [2 a) 
inset, also see Supplementary Material]. Figure |3 a) shows 
the PSDs of the control experiments with a % mask with dif- 
ferent time constants i K =2 ms, 600 /is, 200 /is, 108 /is, and 
72 /is, respectively. All measured PSDs agree very well with 
theory (see Supplemental Material for details). 

We now turn to the measurement of the parity switch- 
ing rate T^Je of the qubit, showing in Fig. |3b) the PSD ob- 
tained without the % mask. The measured spectrum is almost 
flat and white-noise like, suggesting fast qubit parity switch- 
ing. An upper bound can be placed on rff by comparing 
with the theoretical predictions for the observed duty cycle 
Podd — f /(T + Z e ) — 56% and different time constants, T 
and % e . We note that the theoretical predictions in Fig. |5jb) 
assume an RTS amplitude A = F, almost twice as large as 
that in the calibration experiment, because the RTS is just 
from parity switching (see the Supplemental Material). Even 
at x ue = T + T e = 20 /is, there is still some deviation of 
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FIG. 5: (a) Power spectral densities of control experiments using a 
n mask with different time constants z K = f^o + T^i. %o = ^jti m:e 
the time constants associated with the K mask. The color curves are 
data and black curves are theoretical predictions. Inset: schematic 
of the control experiment with a n mask, an RTS with specified time 
constants and generated by an FPGA. The expected RTS amplitude 
is reduced to A con = A K — Aq = P oc idF ', where A n and Aq are the sta- 
tistical averages of the readout during the n mask assuming parity 
switches fast compared to %, P ^ is the probability of the qubit to 
be in the odd parity state, and F is the readout fidelity, (b) Measured 
PSDs and theoretical predictions for the direct experiment on Y% 
without the it mask. Blue curve is the measured Fourier spectrum 
and the color curves are theoretical predictions for different qubit 
parity switching times T oe = T a + T e , where T and T e are the odd- to- 
even and even-to-odd switching time, respectively. The red dashed 
line is the white noise spectrum in the limit T c , T c < t s , the sampling 
time. The near absence of any deviation from white noise at high fre- 
quencies indicates T ,T e < 10 /is. Inset: The PSD at low frequency 
does not show a Lorentzian plateau. 



the theoretical prediction from data. We infer from this that 
qubit parity switching time is faster than our repetition time, 
T ~ T e < 10 /is. This upper bound has been lowered by more 
than three orders of magnitude from the previous estimate in 
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In Fig. HJb), there is a small rise in the spectrum 
below ~ 1 kHz, but its amplitude is an order of magnitude 
smaller than the theoretical prediction at x oe =2 ms with the 
same corner frequency, excluding quasiparticle tunneling as 
the cause. A slow background charge motion can contribute to 
this low-frequency plateau. An analysis of the spectrum down 
to 0.1 Hz with the absence of Lorentzian plateau [Fig. |3b) 
inset] confirms that no slow RTS process is missed. 

The upper bound of z ~ x e < 10 /is places an upper bound 
on the quasiparticle-induced qubit relaxation time 1 /TfLn- 
Recent theory [16] predicts that the quasiparticle-induced 
qubit decay rate is much slower than the even/odd parity 
switching rate. For current device parameters, rff ~ 1/t and 

< 



300 /is. Combining this ratio with the measured T\ = 2 /is 
gives a lower bound T ~ T e > 100 ns. These results do not 
allow us to determine whether quasiparticles dominate the en- 
ergy relaxation of current transmon qubits. However, it is ev- 
ident that quasiparticle tunneling will need to be reduced in 
order to extend T\ beyond ~ 100 /is. 

We have engineered the bandgap profile of transmon qubits 
by combining oxygen-doped Al for tunnel junction electrodes 
and clean Al as quasiparticle traps in an attempt to reduce the 
qubit relaxation due to quasiparticles. However, the measured 
qubit relaxation time is found to be insensitive to the size of 
the superconducting gap and to quasiparticle traps. On the 
other hand, the nondegradation of the qubit relaxation time 
suggests no measurable change of the junction quality from 
oxygen-doped Al. Non-equilibrium quasiparticles leading to 
the charge-parity switching have been observed in all devices. 
Moreover, the quasiparticle-induced parity switching is shown 
to be faster than 10 /is, an upper bound limited by the detec- 
tion bandwidth. This fast parity switching rate, different from 
the results of other experiments in SCPTs, could be due to 
three things. First, the size of the transmon electrodes is much 
larger and the quasiparticle density may depend on the elec- 
trode size. Second, the qubit readout might generate quasi- 
particles or stimulate the tunneling of quasiparticles between 
electrodes. A third possibility is that the interface between the 
oxygen-doped and clean Al layers might not be transparent 
enough. Although our results neither prove nor disprove that 
non-equilibrium quasiparticles limit T\ in transmon qubits, 
it does indicate that quasiparticle-induced energy relaxation 
must be reduced in the future to achieve T\ much longer than 
100 /is. 
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A. Measured Devices 



All transmon qubits presented in the main text are measured by a coplanar waveguide cavity 
in a conventional circuit quantum electrodynamics architecture yj] and the measurements are per- 
formed in a cryogen-free dilution refrigerator with a base temperature of about 20 mK. All devices 
are fabricated on a sapphire (AI2O3) substrate. Optical lithography followed by fluorine-based re- 
active ion etching is used to pattern the coplanar waveguide structure onto sputtered niobium. 
Electron-beam lithography is used to pattern the transmon qubits, which are deposited using a 
double-angle evaporation process. Table U lists the parameters of the transmon qubits discussed in 
the main text. Qubits C, D, E, and F are on the same chip sharing a common readout resonator. 
All qubits have a quality factor Qf ranging between 50,000 and 90,000. 

B. Extraction of Eq, E\, and P ot jd from Histograms 

Due to the switching of the qubit between the even and odd parities, the readout fidelity cannot 
be extracted directly from the measured histogram at one parity peak, but can easily be obtained 
by combining histograms at both parity peaks. For example, Fig. SI shows the cumulative prob- 

TABLE I: Parameters of the transmon qubits. A is the superconducting gap of the electrodes of the Joseph- 
son tunnel junction, f™" x is the maximum qubit transition frequency between ground and the 1st excited 
state, E™ ax is the maximum Josephson energy, E c is the charging energy, g/2n is the coupling strength, f c 
is the bare cavity frequency, Q c is the quality factor of the cavity, and Qf is the quality factor of the qubit. 
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FIG. 1: Histogram of the readout at the even parity peak in Fig. 3(b) of the main text. A threshold 
V t h = 19 mV is chosen to digitize the readout. Blue and cyan curves are the cumulative probability of 
the histograms after digitization. H even — L even , the cumulative probabilities difference below the threshold, 
gives the visibility. The cumulative probability of the histograms at the odd parity peak [Fig. 4(b) of the 
main text] are similar to those at the even parity peak here and are not shown. 

ability after integrating the histograms taken at the even parity peak. H even and L even axe the two 
cumulative probabilities below the readout threshold. Similarly, the counterparts of H dd and L M \d 
at the odd peak are also measured. H even , L even , H dd, and L D dd are related to £o, £i , and P dd (read- 
out false positives, false negatives, and the probability of the qubit at the odd parity, respectively) 
through the following equations: 



L dd = Podd£\ + ( 1 - Podd) (1 - fib) , ( 1 ) 

Leven = (1 -P dd)£\ +Podd(l ~ fib), (2) 
Heven = H od d = \ — Eq. (3) 

Solving these algebraic equations and using the measured values H dd = 0.71, L dd = 0.47, 
Heven = 0.71, and L even = 0.53 gives £q E\ = 0.29 and P dd = 0.56. 



C. Spectrum of an RTS with Finite Readout Fidelity 

In our control experiment, a % mask is used to randomly turn on and off the % pulse applied 
to the qubit. This can imitate the success and failure of % pulse on the qubit and allows us to test 
the sensitivity of our experiment to the qubit parity switching. For simplicity, here we again only 
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consider the case when the n pulse is applied at the odd parity frequency. If the time constants of 
the K mask T^o = i K \ are large compared to qubit parity switching time [which is true even for the 
shortest x k q = 36 /is in Fig. 5(a) based on the measured result x ~ % e < 10 /is from Fig. 5(b)], 
during T^o {l K \) there are enough qubit parity switches to have statistical averages [Fig. 5(a) inset 
of the main text]: 

A K = P odd (l-e 1 ) + (1-Podd)eo, (4) 
A = e . (5) 
Therefore, the n mask will generate an RTS with an amplitude 

A con =A K -A Q = P odd {\ - d - e ) = P odd F < 1 . (6) 

Similarly, for the qubit parity switching rate T$e measurement (Fig. S2), the RTS amplitude is sim- 
ply A = Ai — Aq = (1 — £\ — eo) = F because of the absence of the it mask. Note: the total variance 
of the readout signal in both control and qubit parity switching rate measurement is independent 
of x % and T oe , but depends only on e , £\, and P odd . Each variance is thus constant. Explicitly, 
the total variance of the control experiment is a^ on = [1 — (A K +Aq)/2] (A k + Aq)/2 and the total 
variance of the parity switching rate experiment is Og Xp = { 1 — [P odd M + ( 1 — Podd)Ao] } [P ( ,ddM + 
(1 — P odd )Ao\. Therefore the total areas underneath different curves (the total spectral power) in 
Figs. 5(a) and 5(b) remain conserved, respectively. 

The finite readout fidelity is only expected to generate a white noise background. Therefore, 
the theoretical predictions of the power spectral density (PSD) S of an experiment with an RTS 
would be a sum of the spectra of an RTS and a white noise: 

S = Srts + 2t s of wise , (7) 

where Srts is the PSD of the RTS, t s is the sampling time (10 /is), and <J% oise is the total variance of 
the noise due to the finite readout fidelity. We have used the above equation to get the theoretical 
curves in Figs. 5(a) and 5(b) of the main text. Due to the finite measurement bandwidth, aliasing 
up to three times the Nyquist frequency in the control experiment and up to six times in the qubit 
parity switching rate measurement has been included in the theory curves. However, in the limit 
of a fast RTS signal, x„ ~ x e < t s , the RTS cannot be resolved any more and only a white noise is 
expected with an amplitude S = 2t s o^ xp (again to maintain the total spectral power conserved) as 
indicated by the red dashed line in Fig. 5(b) of the main text. Note: the good agreement between 
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FIG. 2: Schematic of the qubit parity switching rate measurement. It is similar to the inset of Fig. 4(b) in 
the main text except that the RTS signal here is generated from the qubit parity switching. Due to the finite 
readout fidelity, the RTS amplitude is lowered to A = F. 



the theoretical predictions and the control experiment also confirms the assumption that the parity 
switching time is small compared to the time constants of the n mask and the parity switching 
only contributes to the white noise background because of the limited detection bandwidth. 

To get the instrumentation background noise and to ensure that it does not add correlation to the 
measurement, we perform readouts when applying neither the n mask nor a n pulse to the qubit 
and present its spectrum in Fig. S3. In this case, because the qubit is always in the ground state 
and the readout is always performed at /oi ~ 7 GHz where the charge dispersion is negligible, the 
qubit parity switching does not make a difference in the readout and does not generate any RTS. 
Indeed, the spectrum is very close to the expected white noise background Sbg = 2? s £o(l — Co) 
only due to the finite readout fidelity, except for small non-idealities like the \/f noise. 



[1] A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R.-S. Huang, J. Majer, S. Kumar, S. M. Girvin, and 
R. J. Schoelkopf, Nature 431, 162 (2004). 
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FIG. 3: Spectrum of readout with no n pulse applied to the qubit. Blue: measured spectrum in linear scale 
with the qubit always in the ground state; black: expected white noise background due to finite readout 
fidelity. The measured spectrum is very close to the expected white noise background, demonstrating almost 
no correlation in the readout. 
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